The Best Circulant Preconditioners for Hermitian Toeplitz Systems

In this paper, we propose a new family of circulant preconditioners for ill-conditioned Hermitian Toeplitz systems Ax = b. The preconditioners are constructed by con-volving the generating function f of A with the generalized Jackson kernels. For an n-by-n Toeplitz matrix A, the construction of the preconditioners only requires the entries of A and does not require the explicit knowledge of f. When f is a nonnegative continuous function with a zero of order 2p, the condition number of A is known to grow as O(n 2p). We show however that our preconditioner is positive deenite and the spectrum of the preconditioned matrix is uniformly bounded except for at most 2p+1 outliers. Moreover the smallest eigenvalue is uniformly bounded away from zero. Hence the conjugate gradient method, when applied to solving the preconditioned system, converges linearly. The total complexity of solving the system is therefore of O(n log n) operations. In the case when f is positive, we show that the convergence is superlinear. Numerical results are included to illustrate the eeectiveness of our new circulant preconditioners.